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The 0{n) vector model with logarithmic action on a lattice of coordination 3 is related to a 
gas of self-avoiding loops on the lattice. This formulation allows for analytical continuation in n: 
critical behaviour is found in the real interval [—2, 2]. The solution of the model on random planar 
lattices, recovered by random matrices, also involves an analytic continuation in the number n of 
auxiliary matrices. Here we show that, in the two cases n = —1, —2, a combinatorial reformulation 
of the loop gas problem allows to achieve the random matrix solution with no need of this analytical 
continuation. 
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I. INTRODUCTION 



dimensions, so that 



Because of conformal symmetry, universality in 2- 
dimensional statistical systems on a regular lattice is a 
vifell-undcrstood topic To perform a further aver- 
age over the ensemble of random planar lattices, in some 
cases, provides even deeper results, because of the still 
larger symmetry (the discrete version of general covari- 
ance) HHIi. 

The generating function of statistical configurations 
over random lattices can be written as the Feynman 
expansion of a proper action: the replacement of real 
or complex bosonic fields with N x N symmetric- or 
hermitean- matrix fields allows to count lattices of genus g 
with a weight proportional to 7V~^, and a large- A'^ limit, 
achieved via steepest descent or continuous approxima- 
tion of matrix spectra, gives the restriction to planar lat- 
tices (see Q for a recent pedagogical introduction). 

Since the pioneering works on Ising and Potts Models 
0, almost all interesting statistical 2-dimensional prob- 
lems have been studied also in their variant on random 
planar lattices. The specific case of 0{n) vector model 
with logarithmic action on coordination-3 lattices, intro- 
duced by Nienhuis on the regular hexagonal lattice jQ, 
has been solved in the random case by Kostov |^, and 
further studied in particular by Il0| . 

For a given lattice A, with V vertices, the 0(n) parti- 
tion function is here defined as 



2'^(|sf|-l) 



n(i+/3"s.-s,) (1) 



where J7„ is the surface of the unit sphere 5" ^ in n- 
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r(t) 



(2) 



and it is introduced in order to have Z„(A, 0) = 1. 
0(n)-invariance implies 



(.M.f))^A5„,,; (3) 



(4) 



Note that higher-order momenta are not relevant, as they 
do not appear in the polynomial expansion of this 
being a specific feature of the logarithmic action and the 
choice of the coordination number. 

In each monomial of the expansion, each edge 
can contribute either with a factor 1, or with a factor 
/3ns-"^sj"', with a € {!,..., rt}. In this last case, we 
interpret an edge (?,j) as marked, and labeled by a. 

In order to have a non-vanishing contribution after 
integration, each vertex must have an even number of 
marked edges of each species stemming from it, thus, 
as the coordination is 3, the only contributing config- 
urations are the ones of self-avoiding unoriented loops, 
coloured with an index a. Summing over possible loop 
labelings leads to a factor n per loop, and we are left 
with a sum over configurations L of self-avoiding unla- 
beled unoriented loops 



L 



(5) 



where k(L) is the number of loops and \L\ is the number 
of edges in L. In this form of loop gas, the problem 
has a natural analytic extension to complex values of n. 
Note that we can alternatively sum over oriented loop 
configurations L* 



Z„(A,/3) = ^(|) 



fc(L*) 



f3 



|L*| 



(6) 
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Figure 1; A portion of a typical configuration of loop gas on 
a random 3-graph. 



A typical configuration over a random planar lattice looks 
like the one in figure 1. The model shows a nontrivial 
critical behaviour in the range — 2 < n < 2, the cases 
n = 1 and n = 2 corresponding respectively to Ising and 
XY Model. In the case n — > only one loop survives 
and when [3 — oo this model deals with Hamiltonian cy- 
cles ^1^. The case n = —1 can also be seen as a particular 
q — > limit of the g-state Potts Model ■ Recently, it 
has also been connected with the generating function of 
forests on the graph Ts'l. 

In this letter we will show that, in the two cases n — 
— 1, —2, a direct combinatorial reformulation of the loop 
gas problem allows to achieve the random matrix solution 
with no need of analytical continuation in n. 



II. SKETCH OF THE RANDOM MATRIX 
SOLUTION 

Consider the set C of planar graphs A with coordina- 
tion 3. We are interested in the double-series generating 
function 



|A| 



|Aut(A)| 



^n(A,/3), 



(7) 



where |A| is the number of vertices in the lattice. When 
dealing with the 0(n) model, n + 1 matrix fields M, Ei, 
. . . , En are introduced d, 0, 0| . The quadratic part of 
the action is 



Sq = tr 



El 



Ei 



while the interaction part is 



Si = tr 



gP {El 



Ei) M 



(8) 



(9) 



and the partition function (UJ is given by 
Zn{9,l3) = lim -^In / VN(M,{E^})i 



-NiSo+Si) 



(10) 



Notice that auxiliary matrices Ea only appear quadrat- 
ically, although with a non-trivial tensorial form in ma- 
trix indices. At this point one can apply a change of 
variables. Say M = OXO~^, with O G 0{N) and 
A = diag(Ai, . . . , Aat). Change variables from E^ to 
E'^ = 0~^EaO. Integration over M can be replaced 
by integration over A: 



Vn{M) = / d^^AA^(A) 



(11) 



with the Jacobian given by the square Vandermonde de- 
terminant 



A^(A)=ni^'-^^- 



(12) 



We are left with 



Z„(ff,/3)- Jhn^ — In / d^AA2(A) 



2?A'({i^a})e-^('5«+^^), (13) 



with 



n 

'5o = iE-^' + ^EE(^")«^(^- 



(14) 



ij a—l 



i i,j a—l 



We integrate over variables {Ea)ij, and, up to an irrele- 
vant constant, we are left with 

Zn{g,P) = Jii^^^ln / d^AA2(A) 



exp 



i i i,j 

(16) 



This result is, for example, the starting point for all the 
following analysis in Kostov article 0. 



III. ALTERNATIVE COMBINATORICS FOR 
LOOP GAS 

For a given lattice A of coordination 3, call V (respec- 
tively V*) the set of configurations P of unoriented (re- 
spectively oriented) self-avoiding open paths, and |P| the 
number of edges in P. Say that P C L if all the edges in 
P are edges of the loop L, and P* C L* if also orientation 
is preserved. 

Introduce the two partition functions 



Z(A,/?,7) = E/3"^'E7'"-'^' (17a) 

P LDP 

Z*(A,/3,7) = E/3"^'' E ^'"'l^l^'l (17b) 



L'DP* 



Say L (resp. L*) is composed by loops {£fe}, of \ edges, 
then the quantities H17|l are equal respectively to 

Z(A, /3, 7) = E n ((^ + Z^)"*"' - (18a) 

L fc 

Z*(A,/3,7) - En (2(7 + /?)"'=' - 2/31''=') (18b) 

L' k 

that is, each edge of L (resp. L*) can be in P (resp. 
P*) or not, independently, but for the global constraint 
that not all the edges of a given loop are in P (resp. 
P*), which accounts for the subtractions in formulas (|18|l . 
Remark that 

Z(A,/3,-/3) = Z_i(A,/3) (19) 
Z*(A,/3,-/3) =Z_2(A,/3) (20) 
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Figure 2; Contraction of the effective vertex, for the theory 
at n = — 1 (top) and n = —2 (bottom). Thin, thick-solid 
and thick-dashed lines correspond respectively to M-, E- and 
F-edges. 



IV. ALTERNATIVE RANDOM MATRIX 
FORMULATION OF THE PROBLEM 

Now we want to restate the combinatorial theory de- 
scribed in section IfffI averaged over the set of random 
planar 3-graphs, as a random-matrix integral. We start 
from the case n = — 1: 



^-i(.9,/3)-E 



,|A| 



|Aut(A) 



.Z(A, /?,-/?), 



(21) 



Apparently, this theory requires three symmetric matri- 
ces: a matrix M for edges not in L, a matrix E for edges 
in L, but not in P, and a matrix F for edges in P. This 
is not the case: the constraint that F-edges do not make 
loops is not local, and cannot be implemented by a lo- 
cal action. It is necessary to build effective vertices of 
arbitrary degree, containing the connected components 
of P, in order to satisfy the constraint. So, as all P- 
edges are treated explicitly in vertex-combinatorics, they 
do not play a role in Feynman diagrammatic, and both 
the integration and the action only depend over matrices 
M and P, the quadratic part of the action being 



So = tr 



2 ^ 



P" 



(22) 



As in the previous treatment, we have the pure M- and 
P-edge vertices 



5«=tr 



- |Af3 -f g(3E^M 



(23) 



where the factor — /3 in the second term accounts for the 
two external P-legs. This corresponds to vertices of A 
which are not adjacent to edges of P. In the other case, 
say that a path of P has p > 2 vertices (in order to be 
non-empty): the corresponding vertex of the action must 
be built shrinking all the edges stemming from the path 
into a single point. The weight of the effective vertex 
is —{j3gY. It has p external APedges, and two P-edges, 



in a certain cyclic sequence. Summing over all possible 
orderings gives a combinatorics 



E rAEMP'EMP-P' 

p'=0 ' 



(24) 



and thus a further term in the interaction part of the 
action is 



5f =trE 

V=2 



J2 [)jEMP'EMP-P' . (25) 

p'=0 ^ 



Again the change of variables M = OXO~^, E' = 
0~^E0 allows to perform first an angular integration, 
then the gaussian integration of auxiliary variables Py. 
The various terms in the action combine to give 

i i 
1 ^ ^ / \ 



ij p=0 



p'=0 



^E^?"f E^? + ^Er 



(26) 

and a final Gaussian integration gives the quantity l|16|l . 
for n — —1. 

The case n = —2 is very similar. We have to calculate 



Z_2(g,/?) = E 



,|A| 



|Aut(A)| 



Z*{A,(3,-f3). (27) 



Again we have to work only with two matrices, M and P, 
but now, as edges in L are oriented, P is an hermitean- 
matrix field, the quadratic part of the action being 



So = tr 



-M^ + EE 



(28) 
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while the pure M- and E-edge vertices give 



A FINAL REMARK 



(1) 



tr 



-M^ + g(3EEM 



(29) 



For the remaining part of the action, say that an ori- 
ented path of P has p > 2 vertices. The weight of the 
effective vertex is still —{PgY, but now it has p external 
M-edges, one E-edge and one i?-edge, in a certain cyclic 
sequence. Summing over all possible orderings gives a 
slightly different combinatorics 



p'=0 ' 



(30) 



and thus a further term in the interaction part of the 
action is 



Sf^ = tr^(^g)P J] ( ^\eMp' EMP-P' . (31) 

p=2 p'=0 ^ 

Again perform the change of variables, then angular inte- 
gration: the various terms in the action combine to give 



oo P / \ 

+EE(/35)''E !'>'.^f^.'^r^' (32) 

i,j p=0 p'=0 ^ 

= 2 E '^'^ " f E '^^^^ + E ]~ 



and a final Gaussian integration gives the quantity (|16|l . 
for n = —2. 



Calculations of section IIVI can be performed also for 
generic 7 different from —(3. We have 



1 



Zig, /3, 7) - ^Imi^ ^ In / d^A A^X) 



•-p[-^(^E^.-fE^? 

i i 

1 - l3g{X, + Aj) 
-(7 + /3)g(A, + A,) 



(33) 



n I 



with X = 1 for unoriented (resp. x = 2 for oriented) 
loops. More generally, if the weight associated to a loop 
of length £ is — J2f whole procedure can 

be repeated, the last factor in the expression for the ran- 
dom lattice partition function being 



n 



Upi^ - Pf9{X, + A,)) 



(34) 



,ns(l-/3B5(A, + A,)), 

We have denoted the positive weights (3b because we have 
seen how the bosonic matrices £"s naturally give rise to 
the jacobian in the denominator. Weights —(3f could 
be obtained either introducing fermionic matrix fields, 
or by the combinatorial method discussed in this paper. 
Cancellations between numerator and denominator for 
equal values of B are a manifestation of Parisi-Sourlas 
supersymmetry jl4| . 
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